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Abstract. HC1 measurements from HALOE in the northern hemisphere during 
mid-May 1997 revealed vortex fragments in which the chlorine reservoir partitioning 
was strongly pushed toward HC 1 (~ 90 % HC1, ~ 10 % C10N0 2 ), similar to partitioning 
previously observed in the Antarctic vortex region. In contrast, observations of 
C10N0 2 and HC 1 in the northern polar spring, 1992 , and in other years, show these 
species established the balance typical for gas phase photochemical reactions in this 
region (~ 60 % HC1, ~ 40 % C10N0 2 ). Annually, chlorine reservoirs in the winter lower 
stratosphere polar vortex are converted to chlorine radicals via heterogeneous reactions 
on particle surfaces at very cold temperatures (less than about 200 K). As temperatures 
warm in spring, the heterogeneous processes become insignificant compared with gas 
phase reactions, and the chlorine reservoirs are reformed. Measurements through the 
northern winter/spring in 1992 show rapid formation of C10N0 2 , followed by steady 
loss of C!ON0 2 and increasing HC1. Although C10N0 2 measurements are not available 
for 1997 , the HC1 increase in 1997 is observed to be much more rapid and the eventual 
HC1 mixing ratio is about 50 % greater than that of 1992 . 

The observations are examined through comparison with the Goddard three- 
dimentsional chemistry and transport model. This model utilizes winds and temperatures 
from the Goddard Earth Observing System Data Assimilation System and a complete 
integration scheme for stratospheric photochemistry. Analysis of the evolution of HC1 
and C10N0 2 shows that the observed difference in the overall rate of HC1 formation is 
explained by the sensitivity of the gas-phase chemistry to the ozone mixing ratio and the 
temperature. The results show that the model accurately simulates HC1 and C10N0 2 
evolution during these two winters. Model validity is further supported by comparisons 
with O 3 and reactive nitrogen species NO and N0 2 . This analysis provides a sensitive 
test of the lower stratospheric chlorine photochemistry, particularly because the analysis 
considers constituent evolution at a time when the HC1 and C10N0 2 are far from a 
photochemical stationary state. 
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1. Introduction 

The northern hemisphere polar vortex was unusually cold and persistent during 
late winter and spring of 1997 [Coy et al, 1997]. Ozone observed by the total ozone 
mapping spectrometer (TOMS) on the Advanced Earth Observing Satellite (ADEOS) 
was anomalously low [Newman et al., 1997]. Temperatures were cold enough that 
heterogeneous reactions involving chlorine reservoir gases hydrochloric acid (HC1) and 
chlorine nitrate (C10N0 2 ) were expected to dominate the partitioning between chlorine 
reservoirs and radicals. In March 1997, low mixing ratios of HC1 were observed at high 
latitudes by the Halogen Occultation Experiment (HALOE) on the Upper Atmosphere 
Research Satellite (UARS). The low HC1 mixing ratios were inconsistent with gas phase 
photochemistry but consistent with heterogeneous chemical reactions which produce 
chlorine radicals Cl and CIO at the expense of the chlorine reservoirs HC1 and C10N0 2 
[Muller et al . , 1997]. Other articles describing the meteorology and constituent behavior 
for winter/spring 1997 are found in the November 1997 issue of Geophysical Research 
Letters. 

Once the lower stratosphere warms and heterogeneous reactions cease to play a 
significant role in the partitioning among chlorine species, the chlorine radicals are 
converted to reservoirs via gas phase photochemical reactions. One pathway for removal 
of chlorine radicals is combination of CIO with the nitrogen radical N0 2 , In the high 
latitude winter lower stratosphere, most nitrogen radicals are converted to nitric acid 
(HN0 3 ) by heterogeneous reactions, thus the formation of C10N0 2 follows release of 
N0 2 from HN0 3 by photolysis or reaction of HN0 3 with OH. In the northern lower 
stratosphere the increase of C10N0 2 is usually much more rapid than growth of HC1. 
Production of HC1 takes place mainly through reaction of methane (CH 4 ) with atomic 
chlorine (Cl), whose concentration is small enough that this reaction is not favored, and 
the equilibrium between HC1 and C10N0 2 is established slowly. Theoretical studies 
[e.g., Prather and Jaffee , 1990] and analysis of aircraft, ground based and satellite 
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observations from northern winters 1991-92 and 1992-93 support the paradigm of rapid 
loss of CIO and growth of C10N0 2 followed by a slow return to the balance between 
C10N0 2 and HC1 appropriate for gas phase photochemical reactions [see Santee et al , 
1996a, and references therein]. 

The southern hemisphere winter vortex is significantly colder than its northern 
counterpart, and observations from aircraft and satellite indicate the permanent removal 
of HN0 3 from the stratosphere. In the absence of a source of nitrogen radicals, along 
with continued heterogeneous reactions, 0 3 loss proceeds rapidly. Once 0 3 falls below 
about 0.5 ppmv, the mixing ratio of Cl increases, and HC1 formation is rapid. This was 
predicted by Prather and Jaffee [1990], and the rapid growth of HC1 for very low 0 3 was 
seen in HALOE HC1 observations [Douglass et al , 1995]. 

In spring 1997, mid-March HALOE observations at high latitudes showed HC1 
mixing ratios as low as 0.3 ppbv at 500K, indicative of ongoing or recent heterogeneous 
reactions. By mid May, some of the HC1 mixing ratios observed by HALOE at 500 K 
were larger than 2 ppbv. Note that the lowest 0 3 mixing ratios observed by HALOE in 
the vortex, about 1.5 ppmv, were much larger than those which accompany the rapid 
HC1 formation seen in Antarctica, about 0.5 ppmv. 

The Goddard three dimensional chemistry and transport model (CTM), initialized 
in mid-March 1997, produced an increase in lower stratospheric HC1 that is consistent 
with these observations. The increase in the model took place over the course of a 
month. Inspection of the 6 northern springs observed by HALOE suggests that the HC1 
growth and high ratio to inferred Cly at 500K in 1997 is anomalous. The cause for the 
rapid 1997 increase in HC1 is investigated here using comparisons of observations with 
the model fields, and calculations of constituent changes along parcel trajectories using 
the same chemical integration scheme as is used in the 3D model. Observations and the 
simulation are contrasted for the northern high latitude springs of 1997 and 1992. 

There are several goals to this study. The first is to explain the observed behavior 
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of HC1, to identify the factors which control the rise in HC1 once heterogeneous reactions 
directly involving HC1 have ended, and to isolate the factors which are different in 
1997 than in prior years. A second goal is to continue evaluation of the CTM. Direct 
comparisons of the model evolution with that observed for several constituents test the 
ability of the model to produce credible transport during the transition from the winter 
circulation, which is dominated by the polar vortex, to the summer circulation. In 
addition, such comparisons test the model photochemical scheme. In particular, these 
comparisons test the model representation of chlorine photochemistry when the chlorine 
reservoirs are undergoing rapid change driven by gas phase reactions, that is, when the 
reservoirs HC1 and C10N0 2 are in the uncommon situation that their partitioning is far 
from a photochemical stationary state. 

The data used in this study are described in the next section, followed by a 
description of the CTM in section 3. An overall comparison of model and observations 
is given in section 4. The model results for the reformation of the chlorine reservoirs are 
presented in section 5. This is followed by a discussion in section 6 and conclusions in 
section 7. 

2 Data 

The evolution of constituent fields from the CTM will be compared with 
observations from three of the instruments on the UARS satellite. A description of the 
UARS mission is given by Reber et al. [1993]. Each of the instruments of interest here 
is described briefly below. 

2a Microwave Limb Sounder 

The Microwave Limb Sounder (MLS) is described by Barath et al. [1993]. Waters 
et al [1996] provide validation for CIO data obtained using the Version 3 algorithm. An 
improved algorithm has been used to reprocess MLS observations, and the reprocessed 
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version 4 data are used here. Briefly, there were two errors in the Version 3 CIO 
retrieval; an 8% scaling error and an additional high bias. The bias was produed when 
climatological HN0 3 greatly exceeded the atmospheric (gas phase) HN0 3 ; this is the 
case in the high latitude winter lower polar stratosphere when HN0 3 is often in a 
condensed phase [Santee et al., 1996b]. MLS reports observations on alternate UARS 
standard pressure levels, including 100, 46, and 22 hPa. At 46 hPa, the pressure nearest 
the 500 K potential temperature surface, the single profile 1-sigma precision is 0.5 
ppbv; the estimated accuracy is 10%. 

MLS provides temperature profiles which are coincident with the CIO profiles and 
are used for vertical interpolation of the MLS CIO mixing ratios to potential temperature 
surfaces. Although MLS measures temperature at 22 hPa and lower pressures; at higher 
pressures (the region of interest for this work), the temperatures are taken from the 
stratospheric analysis provided by the National Center for Environmental Prediction 
(NCEP), formerly the National Meteorological Center (NMC). 

2b Cryogenic Limb Array Etalon Spectrometer 

A description of the CLAES instrument is given by Roche et al. [1993]. The 
CIONO 2 profiles which are retrieved using algorithm version 7 are described by 
Mergenthaler et al. [1996]. The precision is within 15% between 50 and 10 hPa. The 
maximum systematic error is 28% between 100 and 10 hPa. Profiles retrieved using 
algorithm version 8 are used in this work. The version 7 data set encompasses January 
9, 1992 - May 5, 1993. The version 8 data include this time period and also the earlier 
period October 21, 1991 to January 8, 1992. Data are reported on standard UARS 
pressure levels, which include 100, 68, 46 and 31 hPa, a separation of about 2.7 km. 
Comparisons of profiles shows very little difference between Version 7 and Version 8 
CIONO 2 data sets for winter and spring 1992. 

CLAES temperature profiles are used in the retrieval of the constituent mixing ratio 



profiles as well as for vertical interpolation of the CLAES mixing ratios to potential 
temperature surfaces. The temperature validation is discussed by Gille et al. [1996]. 
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2c Halogen Occultation Experiment 

A description of the HALOE instrument is given by Russell et al [1993]. HALOE 
observes vertical profiles using solar occultation. Approximately 15 sunrise and sunset 
measurements take place daily at each of two near constant latitudes. Profiles for 
the four constituents used here are retrieved using algorithm version 18. The vertical 
resolution of the HALOE observations is approximately 2 km. Validation of the HALOE 
ozone is discussed by Briihl et al. [1996]. The systematic plus random errors are shown 
to be about 18% at 40 hPa; the errors decrease at lower pressure but increase to 30% 
at 100 hPa. Russell et al. [1996] provide validation of HALOE HC1. The total error is 
estimated to be 12% at 5 hPa, and increases to 21% at 50 hPa. The HALOE NO and 
N0 2 measurements are described by Gordley et al. [1996]. There are significant errors 
in both NO and N0 2 below 25 km in spring 1992 due to Pinatubo aerosols. The aerosol 
loading has decreased greatly by the following spring and does not affect retrievals 
thereafter. In the lower stratosphere (below 25 km), HALOE N0 2 shows a low bias up 
to 0.5 ppbv. Between 25 and 40 km, agreement is within 10% of other observations. 
HALOE NO below 25 km, in the absence of volcanic aerosol, compares well with other 
observations when the concentration is large. Between 25 and 40 km, HALOE NO is 
consistently lower than correlative observations; the bias ranges from 0 to 15% at 25 km 
to 15 to 30% at 40 km. 

HALOE provides temperature profiles coincident with constituent profiles which 
are used to interpolate the constituent mixing ratios to potential temperature surfaces. 
Temperatures are retrieved for levels above where aerosols affect the signals (35 km). 
Below this altitude, in the region of interest for this work, the temperature profiles are 


taken from NCEP. 



3. Model Description 

The CTM used here is an improved version of that used in previous analyses of 
observations from MLS, CLAES, and HALOE [Kama et al . , 1995; Douglass et al, 1997] 
and NIMBUS 7 TOMS observations [Douglass et al, 1996]. Winds and temperatures 
used for 1992 simulations are taken from a stratospheric version of the Goddard Earth 
Observing System Data Assimilation System (GEOS DAS) described by Schubert et al. 
[1993]. The horizontal resolution of the assimilation fields in 1992 is 4° latitude by 5° 
longitude, and winds are interpolated to 2° latitude by 2.5° longitude for transport. 

The assimilation system used to generate the meteorological fields for 1997 
incorporates several important improvements. The system is described fully in the 
Algorithm Theoretical Basis Document of the Data Assimilation Office [DAO, 1997]. 
The improvements most relevant to the stratospheric analysis are summarized here. 
First, the computational pole is placed at the equator, thus the longitudinal grid 
spacing decreases as the Coriolis parameter decreases. This coordinate transformation 
reduces noise at high latitudes but has little impact on the tropical fields. Second, the 
top level of the general circulation model used in the assimilation has been raised to 
0.1 hPa, compared with 0.4 hPa in the older system. This change and inclusion of a 
parameterization of gravity wave drag produce more realistic behavior at high latitudes, 
most notably increased descent within the polar vortex. The assimilation fields for 1997 
are produced at 2° latitude by 2.5° longitude horizontal resolution. 

The transport and photochemical contributions to the constituent continuity 
equation in the model are solved sequentially (process splitting). Advection is 
accomplished using the transport scheme developed by Lin and Rood [1996]. This 
scheme maintains sharp gradients and appropriate correlations for long-lived constituents, 
and does not produce unrealistic maxima or minima. For 1992, the assimilation winds 
are interpolated to the model vertical grid. For the 1997 integration, the winds are 
mapped onto the model vertical levels using an integration scheme developed by S. J. 
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Lin (personal communication). This mapping improves vertical transport in the mid to 
upper stratosphere, but has little impact on the lower stratosphere. 

The photochemical scheme includes all gas phase reactions thought to be important 
in the stratosphere; rate constants are taken from DeMore et al [1994]. Reactions on the 
surfaces of stratospheric aerosols are included, with aerosol surface area specified based 
on SAGE observations [Thomason and Poole , 1993] for 1992. For 1997 the surface area 
distribution for a volcanically clean atmosphere is taken from the Scientific Assessment 
of Ozone Depletion: 1991 [Albritton and Watson , 1992]. This distribution agrees 
with aircraft measurements in the lower stratosphere to within 30%. Photolysis rates 
are calculated using temperature dependent cross sections [DeMore et al. , 1994] and 
reduced fluxes that are interpolated from a table-lookup based on the detailed radiative 
transfer calculations from the model of Anderson and Lloyd [1990]. The photolysis rates 
calculated in this way compare favorably with the photolysis benchmark which was 
developed as part of the AEAP program [Stolarski et al , 1995]. This photochemical 
scheme is used in the CTM and also to calculate chemical evolution along trajectories; 
these calculations are initialized using 3D model fields. 

For all simulations, the initialization for long lived species and families follows the 
procedure described by Douglass et al. [1997]. CLAES N 2 O observations from mid 
November 1992 were mapped into a three dimensional grid using potential temperature 
as the vertical coordinate and potential vorticity for the horizontal structure. This 
field is mapped into an initial field, using the potential temperature and the potential 
vorticity fields for the desired day. Ozone is initialized using the same mapping 
procedure with MLS observations nearest to the initial day. For the remaining long-lived 
constituents and families, zonal average fields from the Goddard 2D model [Jackman et 
al, 1996] are mapped onto the 3D grid using N 2 0 as the vertical coordinate. The NO y 
initialization assumes that no widespread denitrification or dehydration has occurred in 
the vortex [Santee et al, 1997]. The fields developed in this way show good agreement 
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with other observations; for example, the total ozone calculated from the initialization 
shows the same structure as TOMS observations. 

The first day of the 1991-92 simulation is November 15, 1991. The simulation 
for 1997 was run in two parts. The first was initialized December 15, 1996, and run 
through March 31, 1997. This integration was made to compare the model CIO fields 
with MLS observations throughout the winter. The model was re-initialized March 18, 
1997 and run through the northern 1998 spring and summer to be used in support of 
the Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) mission. 
The earlier integration agrees fairly well with the POLARIS integration during the 
period of overlap. The O3 fields reveal the same horizontal structure, however, the vortex 
ozone from the Dec. 15, 1997 integration is high compared with observations, which 
impacts the partitioning among the chlorine species. This bias is likely at least partly 
due to a low model rate for 0 3 destruction in the vortex, as reported recently by Becker 
et al., [1998] and Deniel et al. [1998]. Because the O3 must agree with observations to 
produce realistic growth and partitioning of chlorine species, this analysis will focus on 
the re-initialized POLARIS integration. 

4. Overall Comparison of Model and Observations 

This work focuses on the evolution of HC1 in April and May for 1992 and 1997. 
For the model to be useful in this discussion, the model must accurately simulate the 
HALOE observations to the same degree for both integrations. Comparisons of HALOE 
observations of 0 3 , HC1, and sunset NO + N0 2 with model values of 0 3 , HC1, and 
NOx= NO + N0 2 + 2N2O5 are given in Plate 1. Model mixing ratios for the date 
and location of each HALOE observation are compared with the observed mixing ratios 
for pressure surfaces between 68 and 31 hPa for April 15-30 and July 11-13, 1992 and 
between 100 and 31 hPa for April 1-5 and June 6-12 1997. The 100 hPa level is not 
shown for 1992 because the HALOE observations are affected by the high concentration 
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of aerosols following eruption of Mt. Pinatubo. As noted in the validation paper, 
HALOE measurements of NO and N0 2 are accurate in the lower stratosphere when the 
aerosol loading is low and the mixing ratios of NO and N0 2 are high [Gordley et al., 
1996]. Only HALOE observations north of 50°N are considered. 

The 1992 comparisons show that even after six months or more of integration, 
the model fields are in generally good agreement with the HALOE observations. In 
April, the model NO x is up to 50% lower than the observations, which may indicate 
a problem in the model balance between the NO x radicals and HN0 3 for high aerosol 
concentration, or too much aerosol input to the model. In July, when the sunlit period is 
longer and the importance of the heterogeneous reactions on sulfate aerosols is reduced, 
the comparison between the model NO x and HALOE NO + N0 2 is improved though 
model NO x is still low. In April, model 0 3 is within 10% of HALOE observations. In 
July 1992 a high bias in 0 3 (~ 12%) is evident at all three pressures. The model HC1 is 
somewhat smaller than the HALOE HC1, in part due to the general problem that model 
total inorganic chlorine (Cly) in the lower stratosphere is about 20% less than the total 
Cly suggested by observations. Overall, these comparisons show that model fields are 
reasonably faithful to the observations for the nine months of integration. 

The comparisons of model and observations for 1997 are also given in Plate 1. 
Although the model was initialized only a few weeks before, the model and observations 
are in good agreement for early April 1997. Comparisons for late March 1997 (not 
shown) are not as good, suggesting that the agreement between the model and data 
improves as the initialization adjusts to the GEOS DAS wind fields. The model NO x 
agrees well with observations at the three highest pressures. At 31 hPa the model 
NO x is up to 25% lower than the HALOE sunset NO + N0 2 . Low model NO x is 
also seen consistently in comparison to POLARIS ER-2 data (R. S. Gao, personal 
communication, 1998). The model is sampled without regard to local time, and the 
diurnal change in C!ON0 2 is not accounted for, but this is not large enough to account 
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on the left. The model shows remnants of the polar vortex (high HC1, somewhat low 
0 3 ) near 60°N between 0°W and 90°W, between 60°E and 90°E, and near 30°N between 
120°E and 180°E. The model also shows an area between 90°E and 120°E where tropical 
air which has low O 3 and low HC1 has been transported poleward. HALOE observations 
on May 15 are near 60°N, spaced by about 24°, and should show large differences along 
the latitude circle. The model at 62°N is compared with the HALOE observations on 
the right side of Plate 2 ; model features in both HC1 and 0 3 compare well with HALOE 
observations. The good agreement and the higher than expected HC1 mixing ratio 
observed by HALOE and produced by the model in the vortex fragments suggests the 
value of further analysis. This good agreement implies not only realistic photochemical 
changes, but also realistic breakup and transport of the polar vortex fragments, at least 
for the fragments sampled by HALOE. 

5. Evolution of Chlorine Partitioning 

In general, the evolution of chlorine species in the model matches that derived from 
observations for the time period when the chlorine reservoirs are being reformed and 
CIO and C1 2 0 2 are decreasing. This is illustrated for the 500K surface in Plate 3 (a) for 
1992 and Plate 3(b) for 1997. The labeled tick marks are the 15th of each month and 
unlabeled tick marks indicate the 1st of the month. The vortex average is calculated 
each day for model ClOx — CIO + 2C1 2 0 2 , HC1, C10N0 2 and Cly. The fractions 
ClOx/Cly, C10N0 2 /Cly, and HCl/Cly are given in the figure. Also shown on the plot 
are daily vortex sample averages for ratios of MLS daytime CIO, CLAES C10N0 2 , 
and HALOE HC1 to inferred vortex Cly. To plot these observations as fractions for 
comparison with model fractions, the vortex Cly is inferred from the average of the 
HC1, nighttime C10N0 2 , and daytime CIO observations. The rationale is as follows. 

On February 12-17, 1992, the model vortex HC1 is similar to the model average over 
the locations of the HALOE observations, which suggests that the HALOE average is 
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representative of the vortex. After UARS is yawed to view the northern hemisphere on 
February 14, the vortex average of model C10N0 2 is nearly equal to the average of the 
model at the locations of CLAES observations. The vortex average of model C10 x is 
also approximated by the average of model C10 x at the locations of MLS observations. 
Because of the consistency of the model vortex average with the average of model values 
as sampled by HALOE, MLS, and CLAES, we estimate the actual vortex average from 
the individual measurements in mid-February. Vortex average Cl y is estimated to be 2.5 
Ppbv at 500 K based on these observations, compared with about 2 ppbv in the model. 

Comparison of the constituent evolution with observations shows good consistency 
with the model evolution for 1992. Note that the measurement sampling may not be 
representative of the entire vortex, particularly for HALOE, and the model average at 
the location of the observations is also shown in Plate 3 for each constituent. When 
CIO is elevated early in winter, the evolution of MLS daytime CIO follows that of 
model C10 x , the sum CIO + 2 CI 2 O 2 . Much of the difference between model C10 x /Cly 
and MLS ClO/Cly is accounted for by inferring C10 x from MLS CIO, assuming 
photochemical equilibrium between observed CIO and C1 2 0 2 . Model C10N0 2 is about 
a factor of two smaller than CLAES data at this time. The reason for this difference is 
not clear. By early February however, the modeled evolution tracks the observations 
for all the Cly species. Vortex temperatures rose above the heterogeneous activation 
threshold (~ 195 K) on about January 30, 1992 [Newman et al , 1993]. By middle 
February, CLAES and model show C10N0 2 /Cl y » 0.6. C10N0 2 peaks in early March, 
and CIO/ Cly is small and nearly constant for the remainder of the integration. Model 
and HALOE indicate HCl/Cly is « 0.25 in early February. The rate of HC1 increase in 
the model accelerates once ClO/Cly reaches its minimum. The modeled and observed 
HCl/Cly are in good agreement in late March and early April, exhibiting similar growth. 
The vortex comparison ends when the vortex breaks apart and we cannot identify vortex 
fragments either in the observations or in the model. 
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Later in spring, the model and observations indicate little variance in HC1 or 
nighttime C10N0 2 at 500K. The model averages on May 5 for the latitude band 
40°N-50°N for HCl/Cly and C10N0 2 /Cly are indicated on Plate 3(a) by the large 
symbols. These are in good agreement with the fractions taken from averages of HALOE 
observations (about 42 °N) and CLAES nighttime C10N02 in the same latitude band 
on May 5, also indicated by large symbols on the figure. This partitioning is similar to 
that reported by Dessler et al. [1995]. 

For 1997 there are many fewer observations for comparison. There are no 
measurements of C10N0 2 during 1997. The vortex average Cly is estimated from the 
mid-May HALOE HC1 observations and the model C10N0 2 /HC1, and compares well 
with Cly inferred from HALOE CH 4 [ Dessler et al., 1997]. The vortex average model 
ClOx/Cly increases to a maximum of about 0.8 in late January. There is a broad peak 
of about 0.6 throughout February, and a rapid decline in early March. The average of 
the MLS ClO/Cly agree fairly well with the model C10 x /Cl y at the location of the MLS 
observations. C10N0 2 /Cly again peaks at about 0 . 8 , as in 1992, but the maximum 
occurs several weeks later. HCl/Cly rises from about 0.2 in early March to greater than 
0.5 by about March 20 ; both modeled and observed HCl/Cly continue to increase. This 
rise during early March is steeper than is seen in 1992, when modeled HCl/Cly is about 
0.25 on March 1 and reaches 0.35 by about March 20 . The low fraction for HCl/Cly 
on March 1 is consistent with the later occurrence in 1997 of temperatures cold enough 
that heterogeneous reactions which reduce the chlorine reservoirs dominate the chlorine 
partitioning. The area of temperatures cold enough for rapid heterogeneous chlorine 
activation declined sharply in late February, 1997 [Santee et al., 1997]. Intermittent 
small areas of temperature below 195 K, however, persisted through March, which 
explains the relatively higher CIO amounts. In mid March 1997, model HCl/Cly is in 
fair agreement with HALOE HCl/Cly and the model partitioning between HC1 and 
C10N02 is slightly shifted towards HC1 compared to the partitioning inferred for 1992. 
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Throughout April 1997, unlike April 1992, the model HC1 continues to increase rapidly. 
The increase is confirmed by HALOE HC1, which shows a change between late March 
and middle May that is nearly identical to the modeled change. Note that this good 
agreement implies realistic breakup and transport of the polar vortex fragments, at least 
for the fragments sampled by HALOE. 

In summary, the comparison between 1992 and 1997 shows that colder temperatures 
and a more enduring vortex lead to a later recovery of reactive C10 x to HC1 and 
C10N0 2 in 1997. In addition, the recovery pathway in 1997 produces much more HC1 
than CIONO 2 and the the ultimate HCI/CIONO 2 in 1997 greatly exceeds that in 1992. 
We will show that the difference in recovery pathways is primarily the result of lower 0 3 
in 1997 along with lower temperature. The 1997 case is similar to that for the Antarctic 
collar during spring of 1992 [Douglass et aL, 1995]. Note that in spite of the differences 
in chlorine partitioning there does not appear to be a feedback on 0 3 , at least for these 
specific cases. The chlorine-induced 0 3 loss for the two cases is nearly the same after 
March 20, consistent with the low ClOx- 

To account for the springtime differences in the HC1 evolution in 1992 and 1997, 
forward trajectories were calculated for March 28 - April 19 for both 1992 and 1997. The 
initial points for the trajectories are the locations of high latitude HALOE observations 
on March 28. The photochemical evolution along these trajectories was calculated, 
initialized from the 3D model fields, and noting general consistency between the 
initialization from the model and the HALOE observations. 

For 1992, all parcels show the same behavior for HC1 (Figure la). The tendency in 
HC1 is small - the maximum change in HC1 over the 20 day period is less than 0.2 ppbv. 
In 1997, parcels fall into two categories (Figure lb). For parcels in the first category, 
the HC1 mixing ratio is nearly steady, and is similar to the HC1 mixing ratio in the 1992 
parcels. For parcels in the second category, HC1 increases by up to 0.55 ppbv over 20 day 
integration. On average, the HC1 mixing ratio at the end of this integration for parcels 
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in the second category is 0.5 ppbv higher than that for the first category parcels. The 
potential vorticity calculated from the assimilation wind fields is on average 60% higher 
for the trajectories of the parcels in which the HC1 increases, evidence that these parcels 
were part of the polar vortex. At the end of the integration, the results of the trajectory 
calculation for constituents including HC1, C10N0 2 , and 0 3 compare well with the 
model fields at the model location nearest the trajectory endpoint. This is illustrated 
by Figure 1(c), which shows the difference in initial and final HC1 calculated using the 
trajectory model and calculated using the 3D model field at the model location nearest 
the initial and final points of the trajectories. This good agreement supports the use of 
the trajectory model for analysis of the processes contributing to the growth in HC1 and 
leading to a final HCl/Cly in ex-vortex air which is much greater than in surrounding 
air. 

6. Discussion 

The principle gas phase production and loss processes for HC1 in the lower 
stratosphere are 


Cl + CH 4 -> HC1 + CH 3 

(1) 

HC1 + OH Cl + H 2 0 

(2) 


Other production processes such as the reaction of CIO with OH and reaction of Cl 
with H 2 , H 2 0 2 , and H0 2 and other loss processes such as photolysis of HC1 are included 
in the model but play a minor role in the HC1 evolution and are not considered in 
the following discussion. Heterogeneous reactions which take place on aerosols, HC1 + 
C10N0 2 and HC1 + HOC1, are also included in the model but are not significant for 



18 


these trajectories compared to gas phase reactions. 

Neglecting transport and mixing, the continuity equation for the HC1 mixing ratio 
Xhci is written 

dAwri 

= Kci,ch4[C1]X CH4 - K H ci,oh[OH]Xhci (3) 

where X C m is the methane mixing ratio. The 48 hour average production and loss 
for all parcels are compared in Figure 2. The maximum and minimum of the average 
production and loss for all 1992 parcels are shown by the shaded bands. The daily 
maximum average values for each of the 1997 parcels are given individually. Production 
for all parcels in 1992 is about 0.02 ppbv/day. For some of the parcels in 1997 (solid 
lines in Figures 1(b) and 2), production is in a similar range. For the parcels in which 
HC1 is increasing, production is larger, up to an average of about 0.05 ppbv/day. The 
loss is similar in both years for most of the integration (Figure 2b). HC1 is somewhat 
higher in 1997 (Figure 1), but OH is somewhat higher in 1992. (The increase in modeled 
OH is due to Mt. Pinatubo aerosols, which speed the reactions H 2 0 + N 2 0 5 2 

HN0 3 and H 2 0 + BrON0 2 — > HN0 3 + HOBr; photolysis of HN0 3 and HOBr leads 
to higher production of hydrogen radicals.) For all 1992 parcels and some 1997 parcels, 
the loss is approximately equal to production. For parcels with increasing HC1, the loss 
(which is proportional to HC1) increases throughout the integration, and is approaching 
balance with the production at the end of the integration. This is in general agreement 
with Plate 3. In 1992, HC1 is nearly constant by April 1. In 1997, the HC1 levels off 
during mid to late April. When HALOE samples the high HC1 air in mid-May as in 
Plate 2, the HC1 is in approximate photochemical balance. 

It is obvious from Figure 2 that the large difference in production is responsible 
for the observed difference in HC1 evolution. Since W C H 4 is comparable in the model 
in 1992 and 1997, the product of the mixing ratio of Cl and the reaction rate K C i,ch4 
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must be greater for the parcels in which HC1 is increasing. Model output shows that the 
rate K C1)CH4 is larger (due to the lower temperature) and the mixing ratio of Cl is also 
larger. Insight is gained by considering the factors which produce this difference in the 
production of chlorine reservoirs. The production in parts per second can be written 


P HC1 - K oih ,c,Xch 4 x [§5] x [cSoNOd X [ C10N °2] 


(4) 


The fractions Cl/CIO and C10/C10N0 2 are replaced in (4) by relationships which 
are derived from principle production and loss processes: 

[Cl] _ Kcio,o[0] + K C io,no[NO] 


[CIO] 


Kci, 03(03] 


(5) 


[CIO] 


l C10N02 


[C10NQ 2 ] K c10) no2,m[N0 2 ][M] 


(6) 


Neglecting Kcio,o[0] compared with Kcio,no[NO] in the numerator of (5), 
substituting (5) and (6) into (4) and regrouping yields 


P HCl — KcH4,C1^CH4 X 


Kcio,no w Jciono2 _ [NO] 


x 


x 


Kci, 03(03] Kcio,no2,m [N0 2 


xl, 


C10N02 


(7) 


The nitrogen radicals are related through principle production and loss processes 


[NO] Jno 2 

[N0 2 ] K NO) cio[C10] + K NO ,03[0 3 ] 


( 8 ) 


Substituting for N0/N0 2 from (8) into (7) yields 
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p tv" Y KciO,NO JciON02 

*HC1 — J^CH4,Ci2LcH 4 X — rz-^r X 


X 


Kci, O3[03] Kcio,N02,M 

Jn 02 


x Xi 


C10N02 


(9) 


Kno,cio[C10] + KnO,03[03 

The term K NO) cio[C10] dominates the denominator of (8) when CIO is elevated. For 
example, for 0 3 ~ 2.4 ppmv, and the temperature between 200 and 210 K, K no ,cio[C10] 
is larger than K NO , 03 [O3] when [CIO] > 0.3 ppbv. Photolysis of HN0 3 and reaction of 
HN0 3 with OH produces nitrogen radicals, and C10N0 2 growth must precede growth 
of HC1. 


CIO decreases steadily as NOx is produced and C10N0 2 increases, and the reaction 
of NO with CIO becomes less important than the reaction of NO with 0 3 in controlling 
the ratio NO/N0 2 . The reaction of NO with 0 3 becomes dominant once the mixing 
ratio of NOx is of the same order as the mixing ratio of C10 x . There is aerosol 
dependence implicit in this condition - higher aerosol loading is consistent with lower 
mixing ratio for NO x (Plates 1 and 2) and higher CIO [Fahey et al , 1993]. The growth 
rate of HC1 in 1992 increases when ClO/Cly has fallen to about 0.1. At this point in the 
model, the mixing ratio of ClOx is approximately equal to the mixing ratio of NOx- 
Once the CIO mixing ratio has decreased through formation of C10N0 2 so that 
NO + 0 3 is the more important term in controlling the ratio NO/N0 2 , the production 
for a parcel on a constant potential temperature surface may be approximated 


P — KcH4,C1^CH4 x 


Kcio, 


NO 


x 


JciON 02 


X 


JN02 


Kci,03Ao3[M] Kcio,N02,M K NO) o 3 Xo3[M] 


x Xi 


C10N02 


(10) 


This expression depends on the photolysis rates, the methane mixing ratio, the 
inverse square of the ozone mixing ratio, and the temperature through the kinetic rates 
and through the density. It applies to both cases after about March 15. The photolysis 
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rates for CIONO2 and NO2 are nearly equal during 1992 and 1997, The methane mixing 
ratio is also nearly constant. As noted in the Introduction, the polar vortex in 1997 
is cold and long lasting, and the temperature difference between 1992 parcels and the 
1997 parcels in which HC1 grows is about 15K. The ozone mixing ratio difference for 
these same parcels is about 20%, i.e., the ozone mixing ratio is 0.4-0. 6 ppmv lower in 
parcels in which the HC1 grows. The production is also a function of the mixing ratio 
of CIONO2 which is near constant for 1992 but is decreasing throughout the integration 
for 1997 (see Plate 3(b)). As the mixing ratio of C10N0 2 decreases, the production and 
loss come into balance, and the HC1 mixing ratio becomes constant (see Plate 3(a)). 

Equation 10 provides a simple means to evaluate the relative importance of the 
factors which contribute to a difference in production, i.e., the 0 3 mixing ratio, the 
temperature, the CH4 mixing ratio, and CIONO2 mixing ratio, and the photolysis 
rates. The factors which themselves have the largest fractional changes are the 0 3 
mixing ratio and the temperature. The fractional changes in the CH 4 mixing ratio 
and the photolysis rates are small compared to those of the ozone mixing ratio and 
the temperature. Neglecting these smaller contributions, combining the temperature 
dependences for all of the reaction rates, using the low pressure limit for the association 
reaction of CIO and N0 2 , and noting that AM/M = 2.5AT/T on a constant potential 
temperature surface, the expected fractional difference in HC1 production, i.e., AP/P = 
(P1997 - P1992)/ (P 1997 + P 1992) x 2, is given by 


AP — 2A0 3 


Oc 


550 


=-^ ~ -=r + 1.6 -S? 


T 


AT 


T 


( 11 ) 


Here, A0 3 is the difference in 0 3 and AT is the difference in temperature between two 
parcels (1997-1992), and the overbar indicates the average. In Figure 3 we show that 
the expression of the right hand side of (11) is a good approximation for the fractional 
difference in the production for two parcels when the production for each parcel is 
calculated from (10). We compare a parcel from 1997 in which HC1 is increasing with a 



22 


parcel from 1992. For the time period shown, AP/P calculated from (11) and AP/P, 
calculated using the production of each parcel evaluated using (10), differ by less than 
5%. At AP/P = 0.67, the production rate for 1997 is twice that of 1992. The low ozone 
is the most important cause of the more rapid increase in HC1 in 1997 and the higher 
HC1 in steady state. The effect of the temperature difference is far from negligible; on 
average the temperature difference is responsible for 30% of the difference in production 
rate. The high HC1 disappears after the low vortex 0 3 dissipates completely in late May. 

7. Conclusions 

The unusually cold stable polar vortex and the low 0 3 observed in spring 1997 led 
to an observable difference in the reformation of chlorine reservoirs in 1997 compared 
with 1992. The northern vortex was more like the southern winter vortex in that 
it lasted well into spring, but the observed 0 3 was not as low as is observed in the 
Antarctic. Nonetheless, relatively low 0 3 and cold temperatures combine to give a 
faster increase in HC1, and the steady stated-partitioning between HC1 and C10N0 2 is 
pushed towards HC1, resembling the partitioning seen in the Antarctic. The 3D model 
calculations are in good overall agreement with the HALOE observations. 

The breakup of the polar vortex, the transport of fragments, and the photochemical 
change within the vortex fragments produced by the CTM are in overall good agreement 
with the observations available for this time period. There have been other examples of 
lower stratospheric transport on similar temporal and spatial scales which point to the 
applicability of transport calculated with assimilation winds. These comparisons are 
unique in that the photochemical changes produced in HC1 and C10N0 2 are substantial 
throughout the periods of comparison. Furthermore, the model allows interpretation of 
the difference in the growth rate of HC1 in terms of other observed quantities, 0 3 and 
temperature. These comparisons show that the model captures the essential components 
of chlorine photochemistry for an evolving system, a much stronger test than obtained 
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by comparing partitioning for constituents near steady state. 

This work expands the view of reformation of chlorine species. In 1992 and 1993, 
the years for which the CLAES instrument observed CIONO 2 in the northern spring 
vortex, the 0 3 mixing ratio was significantly larger than in 1997. The balance between 
HC1 and C10N0 2 was established slowly; the fraction HCl/Cly exceeded the fraction 
C10N0 2 /Cl y by about 0.2 at steady state. This contrasts with 1997, when low 0 3 and 
cold temperatures combine to produce a balance between the chlorine reservoirs within 
the persistent polar vortex which is pushed strongly towards HC1. In 1997, (observed) 
HCl/Cly ~ 0.8-0. 9, implying that C10N0 2 /Cly ~ 0. 2-0.1. This work emphasizes that 
the rate of HC1 production and the balance between chlorine reservoirs established 
following the cessation of heterogeneous reactions on PSCs is a sensitive function of 0 3 
and temperature. 
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Figure 1. (a) In 1992, the HC1 mixing ratio is constant on all trajectories; (b) In 1997, 
the HC1 mixing ratio is constant on one group of the trajectories (solid lines), but increases 
rapidly for the second group (dashed lines); (c) the difference HCl f - HClj calculated from 
the 3D model values nearest the initial and final points of the trajectory (x axis) agrees well 
with the difference calcualted from the trajectory model (y axis), supporting the use of the 
trajectory model in this analysis. 

Figure 2. (a) The 48 hour average production calculated for some 1997 parcels in Figure 
1(b) (solid lines) is comparable to the production for the 1992 parcels in Figure 1(a) (shaded 
area), but is much larger for other 1997 parcels (dashed lines), (b) The loss is similar for 
both groups of 1997 parcels, until towards the end of the integration when the loss (which is 
proportional to HC1) increases for the parcels which have experienced the growth in HC1. 

Figure 3. The fractional difference in production (AP/P = (Pi 997 - Pi 99 2)/(Pi997 + 
P1992) x 2 where Pi 99 7 and P1992 are calculated from (10), and the contributions to that 
difference from O3 and temperature as given by (11) 


Plate 1 . Comparisons of model values at the location of HALOE observations at OOz on the 
date of the HALOE observation for NOx (top row); O3 (middle row); HC1 (bottom row) for 
April 1992 (leftmost column); July 1992 (left center column; April 1997 (right center column) 
and June 1997 (rightmost column); all HALOE observations north of 50°N are considered. 

Plate 2. The May 15 1997 model O3 on the 500 K surface (top left) and the model HC1 on 
the 500 K surface (bottom left); filled circles locate the HALOE observations. The HALOE 
observations for HC1 (crosses, bottom right) and the difference 0 3 - zonal mean) (crosses, top 
right) show much of the same structure that is found in the model fields (solid lines). 

Plate 3. The evolution of the fractions of Cly for HC1, CIO and C10N0 2 as modeled and 
observed by MLS (CIO), HALOE (HC1) and CLAES (C10N0 2 ) for 1992 (a) and 1997 (b). 
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